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ABSTRACT 

In a stationary, general relativistic, axisymmetric, inviscid and rotational accretion 
flow, described within the Kerr geometric framework, transonicity has been examined 
by setting up the governing equations of the flow as a first-order autonomous dynamical 
system. The consequent linearised analysis of the critical points of the flow leads to 
a comprehensive mathematical prescription for classifying these points, showing that 
the only possibilities are saddle points and centre-type points for all ranges of values 
of the fixed flow parameters. The spin parameter of the black hole influences the 
multitransonic character of the flow, as well as some of its specific critical properties. 
The special case of a flow in the space-time of a non-rotating black hole, characterised 
by the Schwarzschild metric, has also been studied for comparison and the conclusions 
are compatible with what has been seen for the Kerr geometric case. 
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1 INTRODUCTION 

From a mathematical perspective, problems in astrophysical accretion fall under the general class of nonlinear dynamics. 
This occasions no surprise, because accretion, after all, describes the dynamics of a compressible astrop hysical fluid, and the 
fundamental governing equations of such a problem are nonlinear in nature IjLandau fc Lifshitall987l ). One of the general 
issues that is addressed in accretion studies is the physics of the accretor itself, whose gravitational attraction sustains the 
global inflow process. This is especially important if the accretor is a black hole, which by its very definition is never amenable 
to any direct physical observation, and, therefore, its properties can only be known by the gravitational influence it exerts on 
the neighbouring structure of space-time. 

If the accretor is a black hole, the infalling matter has to reach the event horizon on a relativisti c scale of velocity , 
and arguments in favour of this contention have, by now, gained widespread currency jBegelmanl Il978l ; lBrinkmannlll980l ; 
Shapiro fc Teukolsky|[l98i ; IChakrabartilll990h . Given physically sensible boundary conditions, this can only imply that at one 
stage the flow of matter will become transonic, i.e. its flow speed will grow from subsonic values to supersonic values, and in 



doing so, it will at some point ma tch the speed of ac oustic propagation in the fluid. T his transition can happen continuously, 
as in smooth transonic solutions ( Chakrabarti 199p| ). or discontinuously, as in shocks jChakrabartilll989l . [l99ol v In accretion 
studies both possibilities have been subjected to concerted investigation. 

Frequently it happens that the transonic feature is exhibited more than once in the ph ase portrait of stationary solutions, 
i.e. the flow will be multitransonic. This is particularly true in axisymmetric rotational flows l Abramowicz fc Zurekll98ll : Fukue 
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critical solutions in the phas e plane of the flow, i.e . they are associated with critical points — alternatively known as fixed 
points or equilibrium points ( Jordan fc Smith! 19991 ). These solutions may even pass through critical points (as, for instance, 



a flow through a saddle point). In this situation much information about various physical properties of accretion processes 
could be gleaned if these critical points are analysed carefully, which is the central objective of this work. 

In this treatment the nonlinear equations describing the steady, in viscid, rotational, axisy mmetric flow in the Kerr metric, 
have been tailored to form a first-order autonomous dynamical system ( Jordan fc S mith 1999). The critical points of the phase 



trajectories of the flow have been identified first, following which, a linearised study in the neighbourhood of these critical 
points has been carried out. As a consequence of this exercise, a complete and rigorous mathematical classification scheme 
for the nature of the critical points has been derived, and it has been argued that the critical points can admissibly be only 
saddle points and centre-type points for the kind of conserved, axisymmetric and rotational flow under study here. While all 
of these are principally the attributes of the hydrodynamical process itself, the influence of the black hole (the agent external 
to the fluid, but driving its flow nonetheless) has also been noteworthy to the extent that its intrinsic rotational parameter 
affects the character of multitransonicity in general and the properties of an individual critical point in particular. This is an 
important result to have emerged from the study. 

As a special case of what has been done in the Kerr geometric space-time, a similar study has also been carried out for the 
flow in the metric of a non-rotating Schwarzschild black hole. The results are entirely in keeping with the usual expectations 
formed on the basis of the conclusions drawn from the analysis carried out for the Kerr metric. 

Global understanding of the flow topologies will necessitate a full numerical integration of the nonlinear equations of the 
flow. This is never an easy task, but equipped with a clear analytical conception of the local features of the critical points of 
the flow, it becomes possible to qualitatively anticipate the nature of the global behaviour of solutions. In fact, if the system 
has no more than two critical points, then with some knowledge of the p hysical boundary conditions of the flow, absolutely 
correct qualitative predictions can be made about the solution topologies ( Rav fc Bhattachariee! 2002 . 2007 ). 

Particularly with respect to accretion on to black holes, this work follows two earlier works in studying critical properties 
of accretion solutions, with the aid of the mathematical tools of a dynamical systems approach. One work was devoted to 
discussing the behaviour of critical poin ts in a steady, axisymm etric, rotational flow, driven towards the central black hole 
in the pseudo-Schwarzschild framework 1 Chaudhurv et al. 2006t ). The second work dealt with transonici ty and its associate d 
aspects in a fully general relativistic, spherically symmetric flow described by the Schwarzschild metric ijMandal et al.ll2007l ). 
The present work on rotational flows in the Kerr metric may be seen as an extension of these earlier efforts, and it gives 
an example of the relative mathematical ease with which otherwise difficult physical questions in general relativistic fluid 
dynamical problems may be confronted. 



2 THE STATIONARY FLOW IN THE KERR METRIC AND ITS FIXED POINTS 



Neglecting the self-gravity of this general relativistic, axisymmetric, inviscid, stationary, compressible hydro dynamic flow in 
the Kerr geometric space-time, th e conserved equation for the specific flow energy, £ , can be expressed as |Andersonlll989l ; 
Barai et al.ll2004l ; iDas et alj|2006h 



£ — hv t , 

which is actually the relativistic analogue of Bernoulli's equation. The specific enthalpy, h, can be defined as 

, p + e 
h = , 

P 

with e, which contains the rest mass density and the internal energy, being further given by, 
P 

7-1 

The pressure, p, is expressed as a function of the density, p, through an equation of state, 
conditions of constant entropy, 5, the speed of sound is defined as 

2 dp 
c B = -5- ■ 

de s 

All of these establish a connection between the density, p, and the speed of sound, c a , as 



(1) 
(2) 
(3) 

kp 1 , from which, under 

(4) 
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with n, the polytropic index, being defined as n — (7 — 1) . This will subsequently lead to an expression of the specific 
enthalpy in terms of c 2 as 

1 



h 



1 — nc 2 

A further definition gives v t as 



v t = 



f(r) 



(6) 



(7) 



m 

as 



which v is the radial three-velocity of the corotating fluid ( Barai et al. 2004 ; Das et al, 2006). The function /(r) is defined 
Ar 2 A 



A 2 - 4\arA + \ 2 r 2 (4a 2 - r 2 A) : 



with A(r) = r 4 + r 2 a' 2 + 2ra 2 and A(r) 



2r + a 2 . In all of these, the fixed parameters A and a are the sub-Keplerian 



specific angular momentum of the flow and the rotating Kerr parameter, respectively. Following these definitions, a final 
expression for the relativistic Bernoulli equation will be derived as 



£ = 



f(r) 



l-nc 2 y l-v 2 ' ^ 
From the continuity condition, the other governing equation of the flow will be obtained as ( Barai et al. 2004 ; Das et all 



2006 



4nA 1/2 Hp 



m, 
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in which the integration constant, m, is the physical matter flow rat e. The height of the thin disc flow, H(r), under conditions 
of hydrostatic equilibrium in the vertical direction, is expressed as ( Barai et alj 12004 ; iDas et al. 2006 
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Making use of equations (JS| and (|10p in equation @, and eliminating the derivatives of c s with the help of equation 
it becomes possible, under the definition that <?i(r) = Ar 4 and g2(r,v 2 ) — \ 2 v 2 



a 2 Vt + o 2 , to arrive at the relation 
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(ii) 



in which /3 2 = 2(7 + and the primes represent full derivatives with respect to r. 

From the form of equation fTTJ it is easy to appreciate that it is a first-order nonlinear autonomous differential equation, 
whose integration will give the integral solutions in the r — v 2 plane. The critical points of these solutions will be derived by 
the simultaneous vanishing of the right hand side of equation fTTJ and the coefficient of d(v 2 )/dr in the left hand side. This 
will give the two critical point conditions as 
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= 0, 



(12) 



(13) 



respectively, with the subscript "c" indicating the values at the critical points. 

To fix the critical points in terms of the flow parameters, it will be first necessary to make use of both equations (|12|l 
and (|13[) to eliminate c 2 c . Following this, some simple algebraic manipulations, with the help of the definition of g2(r,v 2 ) will 
make it possible to express v 2 entirely as a function of r c , and this will be given by 



f'(r c )gi(r c ) 
/(re)si(re) ' 



(14) 



It will then be possible, with the aid of either equation (|12|l or equation (|13|). to express c 2 c as a function of r c only, and all 
these conditions, substituted in equation JSJ, will deliver the roots of r c in terms of £, A, 7 and a. The critical points will, 
therefore, become fixed in the r — v 2 plane. 
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3 NATURE OF THE FIXED POINTS : AN AUTONOMOUS DYNAMICAL SYSTEM 



Quite frequently for any nonlinear physical system, a linearised analytical study of the properties of the fixed points affords a 
robust platform for carrying out an investigation to understand the global behaviour of integral solutions in the phase portrait. 
This is especially useful in the absence of any well- prescribed and general means of solving nonlinear differential equations, 
which, perforce, have to be solved numerically. 

It has already been shown that the stationary, axisymmetric, rotational flow in the Kerr metric can be reduced to a 
first-order autonomous system, and as such it lends itself easily to a dynamical systems study of its fixed points. To do so, it 
should be necessary to decompose equation into two parametrized equations, given by 

L 
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in which r is a m athematical parameter . Since equations (|15[) are autonomous equations, r does not explicitly appear in their 
right hand sides IjJordan fc Sm ith 1999). 

This kind of parametrization represents the first step towards carrying out a linear stability analysis of the fixed points 
of a nonlinear system, and for the present treatment on disc flows in the Kerr geometric space-time, this will give a complete 
classification scheme for the critical points of the flow. In g eneral fluid dynamics problems — all of which are nonlinear 
problems — this approach is quite common ( Bohr et al.ll 19931) . and in the context of accretion st udies (which, in its essence , 
is the study of a compressible fluid flow), this method has been quite ef fectively adopted before (|Rav fc BhattacharieeJ I 2002: 
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of th e general mathematical aspects of this approach ( Mat sumoto et al. 
19891 ). 



Chaudhurv et a l. 2006; Mandal et al. 2007) . Some earlier works in accretion had also made use 
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Making use of the perturbation scheme, v 2 
continuity condition, 
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it is a straightforward exercise to establish a coupled linear dynamical system in the perturbed quantities 5v 2 and Sr. This is 
given by 
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in which /, gi and g2, all of which are contained in the constant coefficients of the perturbed quantities, are to be read at the 
critical points only. Explicitly written, all the newly appeared constants in the coefficients of equations (|17p are to be given as 
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Using solutions of the kind 5v 2 
equations (|17|l . can be set down as 
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exp(f2r) and Sr ~ exp(rir), the eigenvalues of the stability matrix corresponding to 
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with /, 31 and <?2 to be read once again at the critical points. 

The form of equation (|18fl indicates that the critical points can only be either saddle points ( when Q 2 > 0) or centre -type 
points (when f2 2 < 0), which is just what they should be for a system that is conservative in nature ( Jordan fc Smithll999h . The 



properties of a critical point could be ascertained by making use of the critical point coordinates, (r c ,v 2 ), in equation (|18p . to 
find the corresponding value of Q 2 , and more especially, its sign. This knowledge, along with known and physically meaningful 
boundary conditions of the flow, will give a clear idea of the local behaviour of the integral solutions in the vicinity of the 
critical points. If the system is simple enough, i.e. if it has one or at most two critical points, then a complete qualitative 
impr e ssion of the global behaviour of the solutions can be obtained from this simple analytical exercise (|Rav fc Bhattachariee 
2002, 2007). Occasionally the se situations can actually arise in thin disc flows (relativistic or otherwise) for certain values of 



the re levant flow parameters jDas et al.l 2006p. And t his is almost certainly true for relatively simple spherically symmetric 
flows jRav fc Bhattachariee!l2002l ; iMandal et al.ll2007h . 



4 THE SCHWARZSCHILD LIMIT 



The study carried out so far has been in the Kerr geometric space-time. Although this is a very general case of a relativistic 
astrophysical flow, of no less interest — especially from a theoretical viewpoint — is the case of relativistic axisymmetric 
flows in a spherically symmetric metric. This limit — the Schwarzschild limit — is to be achieved by simply making the Kerr 
parameter vanish (a = 0) in equations (JS|, @ and (|10|) . This will make the black hole a non-rotating accretor, and in the 
space-time described by its gravity, the relativistic Bernoulli equation will be expressed as 



£ = 



r 2 (r - 2) 



1 - net V (1 - v 2 ) [r 3 — A 2 (r — 2)] ' 
while the thickness of the disc, H(r), will go as ( Das et al, 20061 ) 
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The equation of continuity will retain the same form as equation Q even in the Schwarzschild metric, and making use of 
equations ((5JI and (|20[) in equation ([9]), it should be possible to recast the equation of continuity in a suitable form as 
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Following this, it will become quite easy to derive an expression for the gradient of solutions in the r — v plane, and this 
will read as 



i 
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with / 2 being given further by the definition 

2r-3 _ 2r 3 - A 2 r + A 2 
~ r (r - 2) r 4 -A 2 r(r-2)' 

The critical points (labelled by the subscript "c") could be read from the critical conditions delivered by equation 
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Figure 1. Regions of multitransonicity for both accretion (marked A) and wind (marked W) in the parameter space of E and A. The 
labels I and indicate the regions of the lone inner and outer critical points, respectively. The plot has been generated for a = 0.3 and 
7 = 1.33. 
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From the form of the critical points in the foregoing expressions, it is quite evident that with the help of equation (|19p . the 
critical point coordinates can be fixed in terms of £ , A and 7. To understand the behaviour of the critical points, equation (1221) 
would have to be parametrized according to the prescription outlined in Section [3] This will entail writing 

d (v 2 ) = 2/3 2 ^/ 1 (r)-2/ 2 (r) 



dr 
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following which, imposing small first-order perturbations about the critical point coordinates, the eigenvalues of the stability 
matrix of the resulting linearised coupled dynamical system, will be derived as 
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(25) 



in which the arguments of the functions /1 and will be r c , with the primes representing full derivatives with respect to 
r, as usual. With each physically admissible value of r c , the nature of the corresponding critical point will be known from 
equation (|25[) . and just like the flow in the Kerr metric, the critical points will be seen to be only either saddle points or 
centre-type points. They could not have been very different, since the treatment on flows in the Schwarzschild metric is anyway 
a special case of the flow in Kerr space-time. Nevertheless, flows in Schwarzschild geometry merit a separate investigation 
in their own right, and indeed, for spherically symmetric flows, this special case affords a ver y clear pedagogical model to 
understand various accretion-related phenomena, with some surprisingly new features revealed |Mandal et al.ll2007l) . 



5 PARAMETER DEPENDENCE OF MULTITRANSONICITY : GENERAL FEATURES 

In the two foregoing sections, it has been discussed that for the disc flow in the Kerr metric (and for its Schwarzschild limit 
as well), the eigenvalues of the stability matrix associated with each of the critical points can be fixed in terms of the relevant 
parameters of the flow £ , A, 7 and a. This obviously implies that the critical behaviour of the flow can actually be determined 
by these parameters. Multitransonicity is one very important critical feature of the flow and its dependence on the parameters 
£ and A in this Kerr geometric system (for fixed values of 7 and a) has been depicted in Fig.[U 

The region marked by O corresponds to a single outer critical point at large length scales. This occurs for low values of A 
and £. It is not difficult to intuitively grasp the reason for this. Gravity drives the accretion process, and this is manifested by 
the growth of the velocity field. In a rotating flow angular momentum acts against the interest of gravity, and depending on 
the strength of the presence of angular momentum, the velocity field will develop accordingly. Criticality in the flow is achieved 
when its velocity matches the speed of acoustic propagation in the fluid. The resistance raised due to the presence of low 
angular momentum in the rotating fluid can, therefore, be overcome even at large distances (where gravity is comparatively 
weak). Besides this, a low value of £ will imply that the compressible fluid is "cold" and as such it cannot offer much of a 
resistance against gravity with the help of its internal thermal effects. And so once again gravity wins easily, with criticality 
developing at large length scales. All of these features are manifested in the region marked O in Fig. [T] 
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Figure 2. Integral solutions of equation representing region A in Fig.[T] The vertical axis plots M 2 (since it is convenient to measure 
the bulk velocity against the local speed of sound) . Both the axes have been scaled logarithmically for better graphical resolution of the 
plot. The inner saddle point, Si, is connected to itself through a homoclinic path. The solutions have been generated for £ = 1.0001, 
A = 3.15, a = 0.3 and 7 = 1.33. 





Figure 3. Integral solutions of equation JTTJ, representing region W in Fig.[T] The outer saddle point, So, is on a homoclinic path. The 
solutions have been generated for £ = 1.0001, A = 3.6, a = 0.3 and 7 = 1.33. 



On the other hand, when both £ and A are high, their resistive effects can only be overcome with the rotating fluid having 
to fall deeper within the potential well, and, therefore, criticality can be attained only at small length scales, in the vicinity 
of the event horizon of the black hole. This feature has been shown in the region (corresponding to higher values of A and £ ) 
marked I in Fig. [1] All of these arguments can be appreciated much more easily for rotational flows in the non-relativistic and 
Newtonian representation. One might expect that the qualitative character of the physics in this regime, should also carry 
over smoothly to the general relativistic case. At least the parameter space representation in Fig. [T] does n othing to m ake one 
believe otherwise, because a similar pattern is also exhibited for non-relativistic pseudo-Newtonian flows (|Da sl l2002h . 

The multitransonic aspect of the flow has been shown by the wedge-shaped region in Fig. [T] This region has been further 
subdivided into two regions, marked by A (accretion) and W (wind). To gain an understanding of the physical criterion 
behind this subdivision, it should be necessary first to go to equation (O, and then substitute p in it in terms of c s , with the 
help of equation ([5]). This process will lead to the defining of a new parameter for the flow, M = (-yk) n m. This newly defined 
parameter is physically understood to be the entropy accretion rate. 

Now multitransonicity in this Kerr geometric flow implies the existence of three critical points, and under the fundamental 
physical requirement of accretion being a process whereby a flow solution should connect infinity to the event horizon of the 
black hole, the three critical points should be such that there would be two saddle points flanking a centre-type point between 
themselves. Going back to Fig. [1] the criterion to distinguish region A (the accretion region) will be that the entropy accretion 
rate, Mi n , through the inner saddle point must be greater than the corresponding flow rate, Mout, through the outer saddle 
point. The exact reversal of this argument, i.e. A^out > Min, will define the region W (the wind region) in Fig. [l] 

The former situation can be understood very clearly from Fig. [2] in which the integral solutions of equation (|lip have been 
drawn under various boundary conditions. The vertical axis plots M 2 (with M, the Mach number, being defined as M = v/c B ) 
against the radial distance along the horizontal axis. The behaviour of the solution passing through the inner saddle point, 
Si, (whose associated value of M is greater than its corresponding value for the outer saddle point, S Q ) is most interesting. It 
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Figure 4. A comparison of the multitransonic regions in the parameter space of £ and A, for two different values of a, and with 7 = 1.33. 
The dotted curves (coloured blue in the online version) are for a = (the Schwarzschild limit). The continuous curves (coloured red in 
the online version) are for a = 0.3. For the latter case, the shift of the multitransonic region towards lower values of A and higher values 
of £ is quite evident. The area of multitransonicity also increases with the increase in the value of a. 



connects the point Si to itself through a loo p, and hence, to invok e the terminology of plane autonomous dynamical systems, 
this solution is actually a homoclinic path ( Jordan fc Smith 19991 ). 

While Fig. [2] depicts the actual behaviour of the phase trajectories corresponding to the wedge-shaped region marked by 
A in Fig.[TJ the converse behaviour of these trajectories, corresponding to the region W (wind) has been shown in Fig. [3] Here 
the governing principle is that the entropy accretion rate pertaining to the outer saddle point is greater than the one for the 
inner saddle point (i.e. Mout > M\n), and it is evident from the plot that in this situation the solution passing through the 
outer saddle point, So, is a homoclinic path. So a general conclusion that can be drawn is that for a physical flow such as the 
one under study here, multitransonicity will imply the existence of a homoclinic path (which is a flow solution that connects 
a saddle point to itself). The quantitative physical guideline to identify such a solution is to first identify the saddle point 
through which the entropy accretion rate, M, is greater than what it is for any of the other ones. 

For accretion in particular, a long-standing understanding has been tha t the flow ha s to take place at the maximum 
possible rate. One could go back to a pioneering work in this subject, in which iBondil (|1952l ) had conjectured that since there 
would be nothing to prevent the accretion process, it might as well take place at the greatest possible rate, implying that the 
flow would be transonic. While this conjecture was made on the basis of the spherically symmetric flow, which has a single 
critical point (a saddle point) , it would still be relevant for multitransonic disc flows. Once, after starting under suitable outer 
boundary conditions, an inflow solution has passed through the outer saddle point, it will then undergo a transition, so that 
the entropy accretion rate is increased further by the solution passing through the inner saddle point. This transition will 
occur through a standing shock, and this process defines an unambiguous path for the infalling matter to reach the even t 
horizon of the black hole. Various studies have dealt with many questions in this regard l|Chakrabartilll989l : |Pas et alJbood ). 
but these details will be beyond the scope of the present study, which is devoted only to the critical aspects of the flow. 

A further interesting issue related to Figs. [2] and [3] is that in broad qualitative terms, one is evidently a reversed image 
of the other. While the solutions in Fig. [2] are characterised by A4i n > Mout through the saddle points, the defining criterion 
for solutions in Fig. [3] is Mout > Min- It will then be reasonable to suggest that the accreting system will go through 
a state in which M ln = A^ou t- Since integral solutions are allowed to intersect only at the critical points of a dynamical 
system ( Jordan fc Smithlll999l ). this will mean t hat for the condition M in — A4 out , the two saddle points in the phase portrait 
will be connected by two heteroclinic paths only ( Jordan fc Smithll999l ). Drawing the phase solutions in this situation, however, 
will entail a tuning of the flow parameters (and the boundary conditions) with infinite numerical precision, something that 
should be quite impossible in practice. In Fig. [T] the curve separating region A from region W depicts the condition for 
heteroclinicity. As a matter of fact, this curve, as well as the two other curves bounding the multitransonic region in Fig. [T] 
can all be viewed as the loci of various kinds of bifurcation points ( Jordan fc Smith 19991 ) . 

So far the discussion has dwelt on the critical properties of the flow for a fixed value of the Kerr rotating parameter, 
a. It should now be instructive to consider how the variation of a affects the critical properties, because, for the case of a 
rotating black hole, apart from its mass (which defines all length scales in the flow), its spin parameter will also leave its 
imprint on the physics of the accretion process. In Fig. [4] the way in which the spin parameter influences multitransonicity 
has been shown. The dotted curves delineate the region of multitransonicity in the A — £ parameter space for a — (i.e. the 
Schwarzschild limit). The continuous curves indicate the multitransonic region for a = 0.3. It is very obvious that the onset 
of multitransonicity takes place at lower values of A for prograde flows (implied by a > 0). Besides this, the multitransonic 
region is also stretched to higher values of £. All of these lead to the understanding that the Kerr spin parameter (for prograde 
flows at least) favourably affects the multitransonic character of the flow. 
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Figure 5. All positive values of Q 2 indicate that the inner critical point is a saddle point. Its variation in the parameter space of £ and 
A has been shown for a = 0.3 and 7 = 1.33. The accretion region is given by the lightly shaded area (coloured red in the online version), 
while the darker region of the surface plot (coloured blue in the online version) indicates wind. 




Figure 6. The critical point in the middle is a centre-type point, as the negative values of C 2 indicate. The dependence of Q 2 on £ and 
A has been shown for a = 0.3 and 7 = 1.33. The colour scheme remains the same as before. 



6 PROPERTIES OF THE FIXED POINTS IN MULTITRANSONIC FLOWS 

Various general aspects of multitransonicity have been considered in detail by now. To derive some quantitative insight about 
the specific properties of an individual critical point, however, it will be necessary to go further. The first thing to do in 
this regard will be to examine the behaviour of the eigenvalues of the stability matrix associated with each critical point. 
This has to be done by going back to equation (|18p . which gives a dependence of fi 2 on the critical point coordinates. These 
coordinates, in their turn, have a dependence on the parameters £, A, 7 and a. Keeping the last two parameters fixed at 




Figure 7. The outer critical point is again a saddle point, since all values of Q 2 are positive. The dependence of f2 2 on the parameters 
£ and A has been plotted for a = 0.3 and 7 = 1.33. 
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Figure 8. For the inner critical point there is a monotonic growth pattern for the eigenvalues of the stability matrix, O 2 , with respect 
to the Kerr parameter, a. The other relevant flow parameters have been fixed at £ = 1.00035, A = 3.05 and 7 = 1.33. Positive values of 
Q 2 indicate that the innermost critical point is always a saddle point. The growth of Q 2 with increasing a, indicates a strengthening of 
the saddle-like feature of this critical point. 
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Figure 9. Negative values of Q 2 indicate that the middle critical point is always a centre-type point. The variation of the eigenvalues 
shows no monotonic behaviour and there is a minimum near a = 0.2. The other flow parameters are £ = 1.00035, A = 3.05 and 7 = 1.33. 

7 = 1.33 and a = 0.3, the variation of Q, 2 with respect to £ and A, has been plotted in Figs. \5\ [6] and [7] for the inner, the 
middle and the outer critical points, respectively. In all these three surface plots the lightly shaded area (coloured red in the 
online version) represents accretion, while the dark area (coloured blue in the online version) represents wind. It will not be 
difficult to appreciate that all the three surfaces in Figs. 0E] and [JJ will have a two-dimensional projection on the A — £ plot 
given in Fig. [JJ 

The sign of Q 2 indicates the nature of a critical point. When Q, 2 is positive, it will imply the existence of a saddle point. 
And so from Figs. 0and[7] with Q 2 being positive all the time in these two plots, it is very much evident that the inner and 
the outer critical points are saddle points for a multitransonic flow. Which is exactly how it should be to make the whole 
accretion process feasible, because otherwise there will be no open path connecting the event horizon of the black hole and 
the outer boundary of the flow (which, mathematically speaking, will be at infinity) . Having made a note of this qualitative 
similarity between these two critical points, the quantitative differences will also have to be stressed. The first difference is 
that the respective values of Q 2 for either saddle point, differ from those of the other by orders of magnitude. The inner saddle 
point behaves more robustly in this regard. A further difference is that while Q 2 decreases slightly with increasing £ (at a 
fixed value of A) for the inner saddle point (as Fig. [5] shows), the trend is quite the opposite for the outer saddle point. Here, 
for a fixed value of A, there is a growth pattern for fl 2 with increasing £ , as Fig. [JJ indicates. This growth is quite noticeable 
when the value of A is small. For high values of £, however, there is a sharp dip. As opposed to both the critical points in the 
extremeties, the critical point in the middle is always a centre-type point, a fact that is indicated by the negative values of 
f2 2 , given in Fig. [6] Once again, like the outer saddle point, SI 2 increases with increasing £ , at a fixed value of A. 

While Figs. [5] |6] and [JJ indicate a dependence of Q 2 on £ and A for fixed values of a and 7, it will also be necessary to see 
how fi 2 varies with a, having all other parameters fixed. This will reveal how the mathematical properties of the critical points 
are affected by an intrinsic physical property of the black hole. This has been quantitatively represented in Figs. [8l [9l and [lOl 
The two critical points in the extremeties are, of course, saddle points, but once again, in quantitative terms, they behave in 
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Figure 10. The outer critical point is always a saddle point, as the positive values of Q 2 indicate. There is a steady decay in the 
magnitude of f2 2 with increasing a, which implies a gradual weakening of the saddle-like properties. As usual, the other flow parameters 
are £ = 1.00035, A = 3.05 and 7 = 1.33. 



a manner contrary to each other. While, for the inner saddle point, Q 2 grows with increasing values of a (shown in Fig. [SJ), 
there is a steady decline in the magnitude of Q, 2 for the outer saddle point (shown in Fig. 1 10 [> . So, physically speaking, while 
the presence of the Kerr parameter augments the properties of the inner saddle point, at the same time it has an opposite 
effect on the outer saddle point. Both these features are, however, manifested in the case of the centre-type point, for which 
Q 2 decreases initially with increasing a, reaches a minimum value (near a = 0.2), and then starts to increase. These have all 
been shown in Fig. [9] 



7 CONCLUDING REMARKS 

While stationary flows are interesting enough, studying the dynamic aspects of the accretion problem reveals new insight 
related to transonicity, specifically its long-time evolutionary properties and its stability under linearised time-dependent 
perturbations. These are relatively less complicated ventures to undertake for black hole accreti on in the pseudo-Schwa rzschild 
regime, which essentially preserves the simplicity of the Newtonian construct of space and time (jchaudhurv et al"1l2006t ). It has 
been shown for accreting systems, both spherically symmetric and axisymmetric, that generating solutions through a saddle 
point needs infinitely precise fine-tuning of the boundar y condition, but the flow easily attains transonicity if its evolution is 
traced through time (|Rav fc Bhattacharieeil2002l . 120071 ). While it may be proposed that this treatment can be extended to 
proper general relativistic flows, it must also be noted that involving explicit time dependence in a general relativistic fluid 
flow is always a formidable mathematical problem, more so if the flow is compressible and rotating. If, on the other hand, 
this could indeed be achieved successfully, then a whole host of interesting new features would emerge. One issue, related 
particularly to multitransonic flows (which arguably will involve more than one saddle point), is the kind of solution that 
the temporal evolution will select, that is to say which saddle point will the flow finally choose to reach the event horizon 
of the blac k hole, and the physical selection criterion thereof. In this connection other questions like variability and chaotic 
behaviour ( Das et al. 20061 ) in the flow might also conceivably be brought to the fore. 
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